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We examine the attosecond electron recollision dissociation of recently demonstrated experi- 
mentally [H. Niikura et al, Nature (London) 421, 826 (2003)] from a coherent control perspective. 
In this process, a strong laser field incident on D2 ionizes an electron, accelerates the electron in the 
laser field to eV energies, and then drives the electron to recollide with the parent ion, causing 
dissociation. A number of results are demonstrated. First, a full dimensional Strong Field Approxi- 
mation (SFA) model is constructed and shown to be in agreement with the original experiment. This 
is then used to rigorously demonstrate that the experiment is an example of coherent pump-dump 
control. Second, extensions to bichromatic coherent control are proposed by considering dissocia- 
tive recollision of molecules prepared in a coherent superposition of vibrational states. Third, by 
comparing the results to similar scenarios involving field-free attosecond scattering of independently 
prepared and electron wave packets, recollision dissociation is shown to provide an example of 
wave-packet coherent control of reactive scattering. Fourth, this analysis makes clear that it is the 
temporal correlations between the continuum electron and wave packet, and not entanglement, 
that are crucial for the sub-femtosecond probing resolution demonstrated in the experiment. This 
result clarifies some misconceptions regarding the importance of entanglement in the recollision 
probing of Djj". Finally, signatures of entanglement between the recollision electron and the atomic 
fragments, detectable via coincidence measurements, are identified. 



I. INTRODUCTION 

Laser-based quantum control techniques can be loosely 
characterized by two general control ideologies: scenarios 
based in a time-domain perspective, such as pump-dump 
[l| control, and scenarios based in the energy domain 
such as multi-path interference Q and STIRAP Q con- 
trol. The quantum control of reactive scattering also fits 
these paradigms, where examples of multi-path interfer- 
ence (energy-domain) [3] and wave packet (time-domain) 
Q schemes for coherently controlled chemical reactivity 
can be found. In this article, the coherent control of 
attosecond electron recollision dissociation is studied us- 
ing both pump-dump and multi-slit interference schemes, 
and we find that this process offers the first experimental 
demonstration of wave-packet coherent control of reac- 
tive scattering 

The recollision experiment that is the focus of this pa- 
per proceeds as follows. Using femtosecond lasers of in- 
tensities around 10 14 — 10 15 W/cm 2 , one ionizes an atom 
or molecule near a peak of the instantaneous electric field, 
accelerating the liberated electron in the laser field, and 
causing the electron to recollide with the parent ion [f|. 
The ionization/recollision process occurs in less than one 
cycle of the laser field, and repeats every cycle when the 
peak electric field is large enough to cause significant ion- 
ization. For Tksapphire laser systems (800 nm) typical 
of recollision experiments, one optical period is 2.6 fs 
and the recollision of the continuum electron wave packet 
with the core lasts about 500 asec. A number of processes 
can occur upon recollision. The continuum electron can, 
for example, recombine with the ion while emitting its 



excess energy as a burst of XUV radiation [7[, or scat- 
ter elastically thereby takinga sub-femtosecond electron 
diffraction image of the ion @] . Alternatively, as studied 
here, the molecule can undergo recollision-induced disso- 
ciation 

D 2 — ► D+ + c- — ► D + D + + e~, 

where the strong field ionizes D2 in the first step. In the 
second step, the continuum electron excites the bound 
electron from the S g bonding state to the E u antibonding 
state thereby dissociating the molecule. This process has 
been demonstrated experimentally, and used to probe the 
vibrational motion of the nuclear wave packet follow- 
ing the initial ionization of D2 on sub-femtosecond time 
scales [H. 

Here, we consider the dissociative recollision process 
from a coherent control perspective. First, a full di- 
mensional quantum model, based on the Strong Field 
Approximation (SFA), is constructed. The model is 
then validated by successfully simulating the experiment, 
which is shown to be an example of pump-dump control, 
consistent with a perspective long-held by the experimen- 
tal NRC group [101 . Second, the scenario is extended 
to bichromatic coherent control using an initial super- 
position of D2 vibrational states. Considerable control 
is demonstrated, motivating future experimental stud- 
ies. Third, field-free scattering of D^" and e _ wave pack- 
ets is studied in order to connect the observed control 
of dissociative recollision with the recently constructed 
theory of wave-packet coherent control of reactive scat- 
tering Fourth, the role of entanglement, previously 
suggested to be connected to the vibrational probing see- 
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nario 0, [TT|, G3, i s investigated. Specifically, by com- 
paring the vibrational probing scenario to similar sce- 
narios using field-free non-entangled scattering states we 
demonstrate that it is the temporal correlations between 
the scattering wave packets, and not entanglement, that 
allow for probing and control in dissociative recollision. 
Nevertheless, entanglement is still present, and we con- 
clude the paper by identifying its signatures, detectable 
in coincidence measurements. 

For coherent control, these results are of interest be- 
cause they identify new experimentally accessible exam- 
ples of both bichromatic control and wave-packet con- 
trol of reactive scattering. For strong field recollision, 
these results are important because they clarify the role 
of entanglement in recollision-based probing techniques 
0, and motivate new strong field recollision control ex- 
periments. 



II. METHODOLOGY 

A. Strong Field Recollision in SFA 

The SFA, a well-known method in strong field physics 
[lU , is here developed for the dissociative recollision sce- 
nario. We utilize the SFA in the length gauge, and work 
in the single active electron approximation in that only 
the action of the laser on the ionized electron is consid- 
ered, while the bound electron and nuclei do not inter- 
act directly with the strong field. In general, there could 
also be strong laser-induced coupling between the bound- 
electronic states leading to strong field molecular effects 
(e.g. bond-softening, enhanced ionization) that would 
affect the nuclear states. However, the dissociative recol- 
lision experiments being modeled used an angle-limited 
detection scheme to selectively measure D + fragments 
coming predominantly from molecules aligned perpendic- 
ular to the laser field @. Strong field molecular effects 
are minimized for this geometry, hence justifying their 
exclusion from the analysis. Note that the language be- 
low is specific to D2, but the formalism is completely 
general. 

The exact solution for the wave function propagation 
can be written in the form 

|*(t)) = E?(t,to)|*0 (1) 

= — if dt'UXM0^(Oe~ i( *'~ to)5o |*i> 
J t 

+ e - i(i - to)fio |*i) 

where Hq is the field free Hamiltonian, Vi,{t) = E(i) • r 
is the laser- matter interaction for the electric field E(t) 
and the electron position?. The operator U(t,to) is the 
full propagator, defined by 

&U(t,to) = H(t)U(t,t ), U(t ,t )=T, (2) 



where H(t) = H + Vt,(i) is the system+laser Hamilto- 
nian, and / is the identity operator. Note that all equa- 
tions are written in atomic units, h — m e = e = 1. The 
vector potential A(t) is chosen to be 

£ 

A(t) = — - sin(ut)x, (3) 

UJ 

thus giving the electric field 
dA(t) 

E(t) = ^ = £0 cos(wt)x. (4) 

The initial state l^) is assumed to be in an eigenstate 
of the field-free system. The last term on the right hand 
side of Eq. ([T]) describes the evolution of the unperturbed 
component of the initial wave function. It does not con- 
tribute to the dissociative recollision channel and is there- 
fore dropped from subsequent expressions. Expanding 
now the propagator U(t,t') appearing in Eq. (fTj) in a 
manner similar to Eq. JT]), but this time expanding in 
the electron-electron interaction V ee , yields 

|tf(t)) = - f dt' [ dt"U(t,t")Vee(t")U(t",t')x 
Jt Jt' 

VLit'^'-^^,,), (5) 

where another benign term, this time related to direct 
ionization without recollision, was again dropped since it 
also does not contribute to the process under considera- 
tion. 

Before applying the SFA to this exact expression for 
the dissociative recollision, Eq. (JSJ) is first connected to 
the physical picture of ionization/propagation/recollision 
6] introduced above. In particular, consider the inte- 
grand of Eq. ([5]) from right to left. First, the initial 
state l^i) propagates field- free until the time t' , when 
it gets a 'kick' from the laser field This is fol- 

lowed by propagation using the full Hamiltonian for a 
time {t" — t') after which the wave function is affected 
by the electron-electron interaction V ee . The wave func- 
tion then propagates to the observation time t. Thus, 
the partitioning used in expanding the full propagators 
has yielded a propagation sequence that closely resembles 
the physics: The active electron first sits in the ground 
state until being ionized by the laser. This qualitatively 
matches the operators in Eq. ([5]) to the right of, and 
including, Vi(t'). Following ionization, the electron then 
oscillates in the continuum driven by the laser field un- 
til it recollides with the core. These steps qualitatively 
match the effect of the operators V ee (t")U(t" , <'). 

With the formalism firmly connected to the physical 
picture, the SFA can now be applied. The central ap- 
proximation introduced by the SFA is to neglect the in- 
teraction between the ion and the continuum electron in 
the intermediate propagator U(t",t') where the electron 
propagates far from the core before returning to the ion, 
and also in the final propagator U(t,t"). This means 
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that once the active electron has been ionized, it acts as 
a completely free electron oscillating in the laser field. 
The only deviation from this idealization arises due to 
the electron-electron interaction V ee (t"). From the point 
of view of scattering theory, this means that after ioniza- 
tion the continuum electron is effectively being treated 
by a Born-type approximation in the resulting recollision 



event, albeit a Born approximation dressed by the strong 
laser field that drives the electron recollision. 

In order to carry out the SEA, complete basis sets con- 
sisting of the field-free molecular states for the nuclear 
component and plane waves for the active electron are 
inserted between the operators appearing in Eq. ([5]) to 
give 



*(*)> = - E I dE I dt ' I dt " f dk f dk ' f dk " 

, jj Jto Jt> J J J 



(6) 



x (0^> ; k' , | , t') | ; k, t'> (0(f) ; k, i' | (t')e-^ ( *' " f o) | , 

where Ei is the energy of the initial state presently assumed to be an eigenstate of Hq (superposition states will be 

considered below), \<f>a^) are the bound nuclear states on the £ s surface with quantum numbers n = (y, </, m), \4>e\) 
are the continuum states with energy E and and angular momentum J on the £„ surface, the electronic continuum 
(plane wave) states are 



(r|M)=(27r) 



-(3/2) i[k+A(t)]-r 



(7) 



and k is the canonical momentum of the free electron oscillating in the laser field. Within the above outlined 
assumptions, and after projecting onto a particular final state |^ ,;k/,i}, Eq. ^ reduces to 



(8) 



»j;k f ,tm)) = -E f dt ' f dt " I dk 

n Jt Jt' J 
e -(i/2)ti, dr[k / +A(r)] 3 e -iB(t-n(^K k/ y'|V ee (t / 0|^ ) ;k,O 
e -(i/2)ft" ^[k+A(r)] 2 e - l B„( t ''-t^ e - l iJ 3 ( t '- t o) (0 (s) ;k ^/|^^ ) |^ ); 

where E n is the energy of the nuclear state n, and the propagation of the continuum electron was handled using the 
Gordon- Volkov solutions [1J| 



(k',t"\e~*ri' " v ^ dT \k,t') = (k'^'le-t'W'" dr[k+A(r)] 2 | kii ^ 

= e -(V2)/ t f^[k+A(r)] 2 5(k _ k / )) 

which are exact solutions for a free electron oscillating in a laser field. The remaining integrals in Eq. ([ 
using the stationary phase method, keeping the stationary points on the real time axis, yielding 

(<j>M f ,k f ,t\*(t)) = - J2 E ^ (k0 ^ , ^ ,) a c (t c ,4)a p (/ p , n ,^-4 3 V l ( / P,n,4 j) ) 

n j=S,L 



(9) 



are evaluated 



where 



tp(ko,t b ,t c ) 



dr[k f + A(r)] 2 - E(t - t c 



2J t 

(t c ,t b ) = e^l\2i:) 1 ' 



dr[k + A(t)] 2 - E n {t c - t b ) - Ei(t b - t ), 



Fit ffit 

-E(t c )[k f + A(t b )} + E(t b )[A(t b ) A(t c )]^ + I p , n -^ 



(10) 



(11) 



(12) 



2tt 



|tc-*6| 



3/2 



(13) 
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ai{I p ,n,h) = 



/p,n|E(t 6 )r 



1/4 



exp 



1 (2/ p ,n) 3/2 

3 |E(t b )| 



(14) 



the state-to-state ionization potentials are J Pjn = — Ei, the stationary phase points for the time of birth (tb), time 
of collision (t c ), and momentum at birth (ko) are given by 



ko = — A(<b) = — sin(ujtb)x 
to 

Li- 1 * (t c — tb) * sin(ujtb) — cos(wtb) — cos(u>i c ) 

2 



\if sm(ojt c )z 



= \-\ [sin(w4) - sin(wi c )] 2 - D E ,n - I P ,n^ 

A LO Ot c 
I 



(15a) 
(15b) 

(15c) 



and De,ti = E — E n is the state-dependent dissociation 
energy. The derivatives of tb with respect to t c are defined 
through Eq. (flBb ). For each k/, Eqs. (| 1 5|) admit two 
solutions, the long (L) and short (S) trajectories, where 

t ( c S) - 4 5) < ti L) - t[ L) . The sum over j in Eq. {TQ]) 
accounts for these two solutions. 

In principle, Eq. (fT0|) is the final stationary phase 
SFA amplitude for the dissociative recollision process. 
However, some additional simplifications relevant to the 
present work are in order. First, note that while the am- 
plitude Eq. (jlOp correctly captures much of the essential 
physics of recollision problems, such as the qualitative 
dependence of the final yields on laser frequency and in- 
tensity, it typically underestimates absolute yields by 1 
to 2 orders of magnitude. This is due to the neglect the 
Coulomb potential during the ionization step. However, 
methods are available to correct the ionization rates, on 
sub-cycle time scales, by adding the Coulomb potential 
in some approximate way fl5j ]. Second, a c (t c ,tb) in Eq. 
(TIT))) has divergences. These divergences are artificial in 
the sense that the stationary phase method is simply not 
valid at these points. Rather, one should use the uniform 
approximation to treat these regions. When treated cor- 
rectly, a c {t c ,tb) is a slowly-varying prefactor (as a func- 
tion of the final state parameters) to the amplitude Eq. 
(fTOjl and affects little but the absolute yield. Since the 
present study is focused on building a novel perspective 
of this process based upon ideas in coherent control, we 
choose to avoid additional computational complications, 
and simply drop the divergent term, recognize that the 
absolute value is nonquantitative and accept that this 
model correctly captures the dependence of the yields on 



the laser and initial state parameters, but gives incorrect 
absolute yields. Similarly, sub-cycle corrections to the 
ionization rate are not included. Finally, dividing Eq. 
(|15b ) by the maximum value of the vector potential £q jto 
shows that the term containing dtb/dt c is proportional to 
Ip/(2U P ) = 7 2 , where 7 is the Keldysh parameter [la ]. 
In the tunneling ionization regime, 7 is assumed to be a 
small parameter, and thus the dtb/dt c term in Eq. ([T5b ) 
can be neglected. In the calculations present below, this 
term was neglected. However, it was checked that the 
results did not change, at the level of detail considered 
in this work, when this term was used. 

Further simplifications can be made to the ionization 
matrix elements using Eq. I|15b) 

(^;k ,t b \VL(h)\%) = (<j>i 3) ;-A(t b ),t b \VL(tb)m. 

(16) 

For all tb, the final continuum wave is (r|(— A(tb)), tb) = 
(27r)"( 3 / 2 )e i [°]- r = constant [see Eq. 0]. Neglecting also 
angular excitation of the nuclei (rotational time scales are 
much longer than the current time scales of interest), the 
ionization matrix element becomes 

(<t>i 9) ;k ,t b \V L (t b )\*i) cx (^)|#)cos(t b ) (17) 

where \(f>n ) are the nuclear states of the D2 molecule and 
\<t)n^) are the nuclear state of the ion. The overlap 

can be recognized as Franck-Condon factors 
modulating the ionization step. 

With these last considerations taken into account, Eq. 
(dU]) becomes 



frMj^tm)) «-EE ^^^a^^-^ii^) (is) 

n j=S,L 

^ (<f$j;kf\S DR \*i), 

I 



where Sdr, defined by Eq. (|18[) is the approximate dis- sociative recollision scattering operator, and the (now ex- 
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traneous) time-dependence of the electronic plane wave 
basis was dropped. The D + kinetic energy spectrum 
W(E) is found by integrating over the final scattering 
electron and angular momentum states: 



W(E)= Yl 

J=l,3,5. 



dk, 



4"^;^ 15^1**) 



(19) 



The equivalent expression when starting with an initial 
superposition of D2 states Yli C*|^») is 



W(E) 



E 



Ed 



} E,J^ K f\ DB, 



.7=1,3,5.. 

Of interest is also the total yield defined as 



W{E)dE. 



(20) 



(21) 



Interfering pathways that can be used for coherent con- 
trol can be seen in Eqs. (JTSj) and ([20]) . The sum over n 
in Eq. (|18p . that is, the vibrational motion of the nu- 
clei between the moments of ionization and recollision, 
offers a means of pump-dump control, and also provides 
the opportunity to probe the vibrational motion of the 
nuclear wave packet excited to the E u state following 
ionization, the focus of the experiment Q. The pump- 
dump perspective is modeled and discussed in Sec. IIII A[ 
the results of which arc in excellent agreement with the 
experiment, thereby confirming the validity of the SFA 
model. The sum over i in Eq. (|20p provides a means of 
control through preparation of an initial vibrational su- 
perposition state before the strong field is applied. This 
control is similar in spirit to traditional bichromatic con- 
trol schemes, and is explored in Sec. IIII Bt There are also 
interferences arising from the so-called long and short 
trajectories [sum over j in Eq. ([15)1 ] . However, these in- 
terferences lead to rapid oscillations of the final yield as 
a function of k/ , and average to zero once the yields are 
integrated over k/. 

Note that in this formulation the momentum k, re- 
ferred to as the free electron momentum, is in reality the 
relative momentum of the ion and the free electron, and 
that the center-of-mass momentum was assumed to be 
zero. In typical strong field scenarios, this subtlety is ir- 
relevant. However, it is important to stress this point in 
the present study since the field-free scattering formal- 
ism, presented in the following section, will make use of 



both the lab frame coordinates of the electron and ion as 
well as the relative and center-of-mass coordinate. 



B. Field- Free Wave- Packet Scattering 

Consider then scattering in the absence of an external 
field. We work within the S'-matrix formalism uj], [TMl 



S : 



ab 



-i f™ Hdt 



3ab 



%T. 



ab 



(22) 



to calculate the transitions from the initial state |a) = 
|0n 1 ;pi;P<) to final state |b) = \^y,pf,Pf), where 
Pi (p/) is the initial (final) momentum of the scattering 
electron, and (P/) is the initial (final) momentum 
of the ion, both in the laboratory frame. Within the 
Born approximation, the transition matrix elements T a b 
to first order in the electron-electron interaction are given 
by 



T, 



ab 



{2n)6(E h - E a ) 6(K f - K t ) U^y k, 



(23) 



n ' 1 I 1 



with energies 



E b 









2mi 


p} 


p2 

f 

1 — 


2 


2m/ 



Er, 



E. 



(24a) 
(24b) 



Here Kj and k, (j = i, f) are the center-of-mass and 



relative momenta of the e-D^~ system 



K, 



Pi + Pj. 
mj + 1 



(25a) 
(25b) 



and mj = 2m p is the mass of the ion. 

Due to the 5(Kf — Kj) term in Eq. (f2"3"| . which rep- 
resents conservation of total momentum, = = K 
is a constant of motion. The scattering problem is then 
solved for each K individually, after which the yields are 
averaged incoherently over the K distribution. First, the 
projection of the scattered wave function onto the final 
state basis is calculated as 



6 



K4 u) j;k/;K|*/)| 2 = (27T) 1 



x (<ftk/ 



E J dki6(E h - E a ) 



(26) 



Vc, 



With the initial translational momenta of both the electron and ion antiparallel in the laboratory frame, as done in 
the following sections, Eq. |26|) can be reduced to 



\{<j>M;-k f ;K\* f )\ 2 = (2nf 



E(VM 



'> (tt) k, 

} E,J> K f 



(27) 



r 



where 



ho 



2(\k f \y2 + E-E n ). 



(28) 



The partial yields of interest are calculated by integrating 
over the relevant final states 



W K (E) 



E 



J=l,3,5. 



dk, 



(29) 



(30) 



For the case of e + collisions, the electron impact 



W{E) = J W K (E)dK, 



and the total yield is given by Eq. (j21 



matrix elements 



*> (u) k, 
°EJ' K f 



V' 



as previously described (Ref. [5| 



<f>n ; ho ) are computed 



III. CONTROL SCENARIOS 



A. Vibrational Probing and Pump-Dump Control 



spectrum, W D (E D ) = 2W(E/2), using Eq. ((Til). The 
factors of 2 arise from the fact that E is the relative en- 
ergy shared by the D and D + ; since the D + are much 
heavier than the scattered electron, the center-of-mass 
energy of the two D + is negligible and the relative en- 
ergy E is shared equally between the kinetic energies 
of the D + . Here D2 was initially in its ground vibra- 
tional and rotational state, and the field strength used 
was £ = 0.065 au (J = 1.4 x 10 14 W/cm 2 ). This field 
strength is used throughout the paper. The results were 
calculated for a single recollision event. Multiple recol- 
lisions would increase the absolute yield, but would not 
change the relative yields across the spectrum. The four 
panels show the D + energy spectrum for the four differ- 
ent wavelengths (800, 1200, 1530, and 1850 nm) used in 
Ref. 0. The simulated results are in excellent agree- 
ment with the experiment 0], thus validating the SFA 
model. Wd(Ed) clearly depends on the wavelength, and 
this dependence reflects the motion of the nuclear 
wave function on the S g surface between the moments of 
ionization t& and recollision i c , both providing a means 
to probe the vibrational motion as well as offering an 
example of pump-dump control. 



Consider first the original experiments of Niikura et 
dl. [H, where the D+ kinetic energy spectrum was mea- 
sured as a function of laser wavelength. The basic idea 
was to exploit the frequency-dependence of the time de- 
lay between the moment of ionization tf, and moment of 
collision t c , implicit in Eq. (|15b ). in order to probe the vi- 
brational motion of the vibrational wave packet. This 
vibrational probing scenario can readily be identified as a 
pump-dump control scenario: the ionization step pumps 
the nuclei to the S ff surface while the recollision dumps 
the nuclear wave packet to the dissociative surface. Vary- 
ing the laser frequency controls the pump-dump time de- 
lay. In fact, the vibrational probing experiments were 
guided by this close analogy with pump-dump control. 

Figure [JJ demonstrates the vibrational probing/pump- 
dump scenario using the SFA model. Here, and through 
the rest of the paper, we plot the D + kinetic energy Ed 



B. Bichromatic Coherent Control 

The orthodox bichromatic coherent control scenario 
utilizes CW lasers and involves first creating an initial su- 
perposition state and then driving the superposition to 
a final state where the components of the superposition 
can interfere. This scenario can be extended to the time 
domain, and to attosecond dissociative recollision, in par- 
ticular by creating an initial vibrational superposition in 
D2. Since the recolliding electron has a broad bandwidth 
(energies of recollision vary from to 3.17 U p where 
Up = (£o/oj) 2 /4: is the pondermotive potential), the var- 
ious initially-populated vibrational states will overlap in 
the nuclear continuum on the antibonding surface follow- 
ing collisional excitation, and make possible the bichro- 
matic control of the D + energy spectrum. 
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FIG. 1: Pump-dump control of the dissociative recollision D + 
spectrum. The ionization event provides the 'pump', while 
recollision provides the 'dump'. The pump-dump time delay 
is controlled by varying the driving frequency. 



This scenario is illustrated in Fig. [2] using a superposi- 
tion of the \v) = |0) and |1) vibrational states of D2, with 
rotation in the ground state. The results correspond to 
a single recollision event and a driving field of frequency 
lo = 0.0569 au (800 nm). Panel (a) plots the final D+ 
spectra starting from |0) (solid) and |1) (dashed) indi- 
vidually, while panel (b) plots the spectra arising from 
the vibrational superpositions |+) = |0) + |1) (solid) and 
|— ) = |0) — |1) (dashed). Panel (c) plots the total in- 
tegrated D + yield Wt for the vibrational superposition 
state \<j)) = |0) + e^|l) as <j> is varied. The range of con- 
trol is significant. These results should strongly motivate 
an experimental study of bichromatic coherent control 
via strong field dissociative recollision: the D+ spectrum 
and total yield can be actively controlled by using phase- 
coherent initial states. 

As in the previous section, these results include only 
a single recollision event, thereby effectively simulating 
the result of driving the recollision with a single-cycle 
pulse. In the case of pump-dump control, adding more 
cycles did not change the resulting D + spectrum. How- 
ever, this is not the case for the bichromatic scenario. If 
the number of cycles is large enough so that the total 
time duration of the driving pulse is larger than the vi- 




E D (eV) 

FIG. 2: Single-cycle bichromatic control of dissociative recol- 
lision. D + spectrum for initial (a) |0) vibrational state (solid), 
and initial |1) state (dashed); (b) initial |0) + |1) superposi- 
tion (solid) and |0) — |1) superposition (dashed), (c) D + yield 
from |0) + exp(i<fi)\l) state as a function of cj>. 

brational period of the initial superposition state, then 
control will disappear, since each cycle would contribute 
to the final D + spectrum with a different phase 4> due 
to the evolution of the vibrational eigenstates from one 
cycle to another. However, control can be restored in 
the many cycle regime by using pulses at two different 
frequencies, 



E 2 (t) = (£ /2)[cos(wii) +cos{u 2 t)]x 
= £q cos((Aw/2)t) cos(wt)x, 



(31) 



where lo\ = lo + Aw/2, lo 2 = lo — Aw/2, and Alo is set 
equal to the energy difference between the two states of 
the initial superposition. In this case, the two frequen- 
cies will cause beats in the envelope of the long pulse 
that are timed to the motion of the internal state super- 
position (see Figs. [3^ and b). Here the central frequency 
is set to u> = 0.0569 au (800 nm). Varying the relative 
phase of the superposition, or alternatively varying the 
relative phase of the two frequencies, then varies the rel- 
ative timing between the internal state motion and the 
envelope beats, and allows one to control which phases of 
the vibrational superposition contribute to the recollision 
events. 

Figure [3J: shows the control, as a function cf>, when 
using a monochromatic field (dashed) and the two-color 
field [Eq. flni>] with lo = 0.0596 au and Aw = E v=1 - 
E v —q = 0.0135 au The calculations included eleven cycles 
of the 800 nm carrier, corresponding to two vibrational 
periods of the D 2 superposition. When solving the sta- 
tionary phase equations for the two-color laser field, the 
field and vector potential were treated as pure sine waves 
over the half cycle, 



E 2 (i) ~ £ cos((Aw/2)r„)cos(a;i)i 



(32) 
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FIG. 3: Many-cycle bichromatic control of dissociative rec- 
ollision. (a) Density of D2 vibrational wave function as a 
function of time. The undulation of the density is highlighted 
by the superimposed wave, (b) Incident laser intensity as a 
function of time for the two-color field of Eq. (|31[) with So =, 
u> = 0.0569 au, and Aco = 0.0135 au (c) Control of D+ yield 
as a function of <j> f° r a monochromatic incident field of fre- 
quency uj (dashed) and the two-color field (solid). 



and 



A 2 (t) « — - cos((Aw/2)T„) $a\(ut)x, (33) 



where the amplitude of the beat envelope cos((Aw/2)T„) 
was evaluated at the time T n corresponding to the peak 
field strength of the half cycle, defined by cos(wT„) = 1. 
For ease of comparison, both control plots were normal- 
ized so that the average yield as a function of <p is unity 
to compensate for the different ionization yield that de- 
pends exponentially on the instantaneous electric field 
strength [Eq. (|14[) ]. As can be seen from Fig. [3fc, the 
control with the monochromatic wave is greatly dimin- 
ished for many cycles, and the control eventually reaches 
zero in the limit of infinite number of cycles. However, 
control in the case of the two-color field persists, and does 
not diminish further when additional cycles are added. 



C. Coherently Controlled Reactive Scattering 

An interesting question arises when one considers dis- 
sociative recollision from the perspective of coherently 
controlled reactive scattering. Traditional (non-control) 
scattering scenarios [13, EH are almost all conducted in 
a field-free environment and use nearly monoenergetic 
beams. However, the present scenario necessarily in- 
volves the presence of a strong laser field and wave pack- 
ets of translational momentum. Coherent control of re- 
active scattering in the continuous beam regime Q is 
now well understood, and the wave packet extension has 



recently been explored Q ■ Although dissociative recolli- 
sion clearly involves a scattering process, it is also cou- 
pled to the ionization step and the in-field propagation, 
and hence it is not entirely clear to what extent it can 
be considered as a true example of coherently controlled 
reactive scattering. The easiest way to address this issue 
is to compare the strong field recollision scenario to its 
D2 + e~ field-free counterpart using the field-free scat- 
tering formalism presented in Sec lIIBl 

The incident scattering state used to mimic the recol- 
lision event has the form 

l*o} =J2 I rfKdka n ^(k,K) \^\ k,K), (34) 

n J 

where a n and iftQc, K) are the vibrational and transla- 
tional expansion coefficients respectively. In the recol- 
lision scenario, the collision between the electron and 
the ion lasts about 0.5 fs. For the laser intensity used 
above, at 800 nm, the peak recollision momentum is 
Pmax = y2 x 3.171/p f=a 1.5 au. This recollision wave 
packet can then be adequately represented in the field- 
free case by the following translational superposition 



V.(k,K) 



1 



V^AKAk 



x exp 



exp 



k-k 
Ak 



K-K 
AK 



(35) 



with relative coordinate parameters given by Ak = 0.5 
au and ko = 1.5 au, parameters chosen to mimic the 
recollision wave packet observed in the simulations of the 
previous sections. Since the scattering is independent of 
the center-of-mass parameters, any values of AK and Ko 
can be used. Here we set these parameters to AK =1.0 
au and Ko = au, where ko and Ko are parallel to 
each other. In order to match the vibrational wave 
packet created in the strong field case, the vibrational 
coefficients are set to 



(f)|4 l) )exp 



3 |E(f fe )| ltj " T \ 



(36) 



where ((j)^ \4>o^) are the Franck-Condon overlaps, and 
the exponential mimics both the ionization amplitude 
[see Eq. (fT4")l ] and the time delay t<j = 2.0 fs between 
the momentum of ionization and recollision (r^ is ^3/4 of 
the laser cycle @). This distribution of vibrational states 
on the surface follows from the above presented SEA 
theory by evaluating Eq. (fT8|) after ionization but before 
recollision, and has been observed experimentally [lj|. 
Crucial from the perspective of coherent control is that 
the combination of a spread in population of vibrational 
states |^n 9 ^) [Eq. ([5rj]) ] and spread in incident momen- 
tum [Eq. ([35]) ] leads to multiple energetically degenerate 
pathways to product states. 

Figure S] plots the resulting D + yields for the field- free 
scattering. Upon comparing panels (a) and (b) of Figs. 
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FIG. 4: Field-free wave-packet coherent control of electron 
impact dissociation. D + spectrum for initial (a) |0) vibra- 
tional state (solid), and initial |1) state (dashed); (b) ini- 
tial jO) + |1) superposition (solid) and |0) — |1) superposition 
(dashed), (c) D + yield from |O)+exp(i0)|l) state as a function 
of 6. 
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S] and [H one sees that the D + spectra resulting from 
the strong field and field-free scenarios are essentially the 
same. Small variations are present, for example the dis- 
tributions in the field-free case arc all a bit narrower than 
their field-free counterparts, but this is likely due to small 
differences between the actual recollision state and that 
used to mimic this state in the field-free scenario. When 
comparing the phase dependence of the D + yields [panel 
(c) of Figs. H] and [2j, one sees that the general sinu- 
soidal dependence is perfectly captured by the field-free 
scenario. It is therefore clear that the strong field and 
field-free cases behave in essentially similar ways, and 
hence that the dissociative recollision scenario offers an 
example of wave-packet control of reactive scattering. 



D. Entanglement in Dissociative Recollision 

We now address the importance of entanglement in dis- 
sociative recollision. Specifically, the literature suggests 
a connection between the sub-femtosccond time resolu- 
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FIG. 5: Non-entangled wave packet probing of vibrational 
motion. To be compared with Fig. [1] The Td values corre- 
spond to the pump-dump time delays for the wave lengths 
used in Fig[T] 



entanglement is not necessary when working with spe- 
cific: scenarios utilizing wave packets of translational mo- 
tion [5j involving energetically degenerate pathways. The 
dissociative recollision scenario features such pathways, 
as noted above, and hence it is conceivable that entan- 
glement is neither crucial for control nor for vibrational 
probing. This is examined below. 

To gain insight into the role of entanglement, we con- 
sider once again the analogous field-free case with inde- 
pendently prepared wave packets of and e~, but this 
time without entanglement present. That is, our initial 
state is 



a connection Derwecn tne suD-iemcosccona time resoiu- r 
tion present in the vibrational probing of the wave 

— 4- lr] h il I -i oil t^i — — — ,i~~~~~4- „ J 



packet and entanglement I 111 . [r3 | . The entanglement 
under consideration is that between the electron and ion 
after ionization and before recollision: when the D2 is 
ionized, the resultant wave function is separable in rela- 
tive and center-of-mass coordinates, similar to Eq. (|34[) . 
but not separable in the ion and electron coordinates. 
The question of entanglement is also central to the co- 
herent control of reactive scattering 0, Q. In that case, 
initial state entanglement is required for control when 
using superpositions of a few isolated translational plane 
waves |4J. Recent results show, however, that initial state 



(37) 



where 

V( P ,P) 



I 

x exp 



exp 



1 / P-Po 

2 V AP 



1 ( P Po 

2 ^ Ap 



(38) 



the electronic parameters are given by Ap = 0.5 au and 
Po = 1.5 au, and the ionic parameters are AP = 1.0 au 
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FIG. 6: Signatures of entanglement in the multiparticle distri- 
butions, (a) Electron momentum spectrum, for Ed ~ 6 eV, 
showing ATI structure, (b) Coincident D + /electron spec- 
trum, for kf yV — 0, showing ATI-like structures in the D + 
kinetic energy distribution. Both spectra are taken at kf, z —Q. 



and Po = 0. The a n are the same as those used above. 
Recall that the Po and po are momenta in the laboratory 
frame. 

Figure[5]shows the results of mimicking the dissociative 
recollision probing/pump-dump scenario using the field- 
free non-entangled scattering scenario. The 'pump-dump 
time delay' was varied by controlling r^, and results for 
the four values of corresponding to the wave lengths 
used in the strong field pump-dump scenario are shown. 
Upon comparing Fig. [5] with Fig. [TJ one sees that the 
non-entangled field-free scattering scenario does capture 
the vibrational motion just as well as did the strong field 
recollision scenario. 

Note that the Wd(Ed) peak height increases as Td be- 
comes larger in the field-free scenario, but the analogous 
peak height in the strong field case remains roughly con- 
stant. The change in total yield in the field-free case 
is present because the cross section for electron impact 
dissociation of depends strongly on the internuclear 
distance; as the bond is allowed to stretch, it becomes 
easier to dissociate. In the strong field case, as the wave 
length is varied and the pump-dump time delay increases, 
the continuum electron wave packet spreads more in the 
transverse direction, thus counteracting the expected in- 
crease in the total yield as the bond length stretches. 

These results directly imply that it is the temporal cor- 
relations between the incident electron wave packet and 
the vibrational motion of the ion that is necessary for 
probing/control, but entanglement is not. Further, the 
general agreement between Fig. [5] and Fig. [1] strength- 
ens the argument that dissociative recollision scenarios 
can be considered as examples of wave-packet coherent 
control of reactive scattering. The theory of wave-packet 
coherent control of reactive scattering is discussed in de- 
tail in Ref.[|. 

Although electron-ion entanglement is not central to 
probing or control, the fact remains that entanglement is 
present. The entangled nature of the multiparticle den- 
sity can be revealed through coincidence measurements, 
where the D + kinetic energy spectrum is measured to- 
gether with the recollision electron. The ionization, and 



subsequent recollision, certainly leads to multiparticle 
correlations in the outgoing state. However, to demon- 
strate entanglement, one must show non-separable wave 
function characteristics, which include correlated observ- 
ables and coherence. When using a driving pulse with 
many cycles, the outgoing flux originating from each half 
cycle will interfere, hence providing a means to verify the 
coherence of the final wave function. 

Figure [6^ shows a slice through the final momentum 
distribution of the scattered electron taken at Ed = 6 
eV and kf :Z = 0, calculated for a five-cycle 800 nm pulse. 
The ring structure, similar to the well known Above- 
Threshold Ionization (ATI) peaks, appear as result of the 
periodic nature of strong field ionization. Each half-cycle 
an electron wave packet is launched into the continuum. 
All these wave packets then interfere to produce the rings. 
In an energy representation each ring is separated by one 
unit of the photon energy. Since dissociative recollision 
also occurs each half cycle, similar ATI-like structures 
exist in the D + energy spectrum, FigJSJa, again a result 
of interferences between the outgoing D + flux originat- 
ing from each half-cycle. The spectrum shown in Fig. [BJd 
corresponds to kf^ y — and kf tZ = 0. Unlike one-electron 
ATI, the ring structures resulting from dissociative rec- 
ollision only appear when both the recollision electron 
and the D + are measured in coincidence, that is, the 
structure exist only in such spectra. Further, since the 
rings are a result of quantum interference, the result- 
ing wave function is necessarily coherent. Observation 
of these ATI-like rings in the coincidence spectra would 
then demonstrate the entanglement that is necessarily 
present during the dissociative recollision process. 

IV. SUMMARY 

We have considered the attosecond dissociative elec- 
tron recollision in D2 from the perspective of coher- 
ent control. Both pump-dump and bichromatic control 
schemes were demonstrated in this system. Further, by 
direct comparison of the strong field recollision scenario 
with analogous field-free cases, it was shown that this sys- 
tem offers an example of wave-packet coherent control of 
reactive scattering. By constructing and analyzing simi- 
lar scenarios involving scattering of non-entangled field- 
free wave packets, it was shown that the electron-ion en- 
tanglement effectively plays no role in the control sce- 
nario presented herein nor in the closely related problem 
of recollision-based sub-femtosecond vibrational probing 
Q. Finally, detectable signatures of electron- ion entan- 
glement in the outgoing multiparticle states, arising from 
temporal interferences, were identified. 
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